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Abstract We conducted molecular phylogenetic analyses of Japanese Rhacophorus species, especially of R. schlegelii 
and R. arboreus from the mainland, based on samples encompassing their known distribution ranges, and discussed 
about evolutionary history of Rhacophorus species within Japan. The common ancestor of Japanese Rhacophorus, 
except for R. owstoni from southern Ryukyus, was estimated to have diverged from a lineage occurring mainly in 
China about 7 MYBP. Both R. schlegelii and R. arboreus are genetically largely divergent between regions of eastern 
and western Japan, and this seems to have been promoted mainly by retreat to refugia. Retreats of the two species to 
different refugia sometimes in the past seem to have led restricted distribution of R. schlegelii in eastern and R. arboreus 
in western Japan, and brought their intraspecific variation patterns in morphology and breeding habit. 
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1. Introduction 


Some amphibians like the newt genus Cynops and the 
frog family Rhacophoridae have Japanese mainland as 
their northern limit of distribution. It is suggested that 
some species of such groups contain populations that 
were affected by climate changes in the Quaternary 
glacier and interglacial periods (e.g., Cynops pyrrhogaster 
[Tominaga et al., 2013]; Buergeria buergeri [Nishizawa 
et al., 2013]). A rhacophorid genus Rhacophorus occurs 
from Southeast to East Asia and contains about 92 species 
(Frost, 2019). In Japan, R. owstoni, R. viridis, and R. 
amamiensis occur in the Ryukyu Archipelago, while two 
species, R. schlegelii and R. arboreus are distributed 
on the mainland, and are northernmost species of the 
genus and family (Matsui and Maeda, 2018). This leads 
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us to expect the presence of genetic variations in these 
species induced by climatic oscillations like Cynops and 
Buergeria shown above. 

Rhacophorus schlegelii is distributed in Honshu, 
Shikoku, and Kyushu Islands, and is characterized by 
breeding under the ground. Unlike R. arboreus, with 
which the species is sympatric in Honshu, it has small 
body size and yellow iris, and usually lacks dorsal 
marking. Rhacophorus arboreus occurs only in Honshu, 
and breeds on trees or on the ground but never in the 
soil. Compared with R. schlegelii, it has a large body 
size, sometimes with dorsal marking, and has reddish iris 
(Maeda and Matsui, 1999). 

A few previous studies reported geographic variation of 
these two species. Okada and Kawano (1924) recognized 
two varieties of Polypedates schlegelii Ginther, 1858: 
P. schlegelii var. arborea with a dorsal marking and P. 
schlegelii var. intermedia without such a marking from 
northern Honshu. Later, Okada (1930) elevated the two 
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varieties to subspecies of R. schlegelii. Now R. schlegelii 
intermedia is synonymized with R. arboreus, which is 
treated as a species distinct from R. schlegelii (Nakamura 
and Uéno, 1963). More recently, Iizuka et al. (1990) 
reported geographic variation of nucleolar organizing 
regions in both R. schlegelii and R. arboreus. Using RFLP 
method, Wilkinson et al. (1996) analyzed 12S rRNA-16S 
rRNA of R. arboreus and reported genetic differentiation 
between eastern and western populations. However, no 
comprehensive studies have been made in either species 
encompassing whole distributional areas or estimating 
formation histories of populations. 

It is expected that elucidation of historical formations 
of geographic variations in these two species would 
help understanding evolutions of their morphology and 
breeding habits. Additionally, it is biogeographically 
interesting that R. schlegelii occurs on Honshu, Shikoku, 
and Kyushu like many other Japanese frogs, but R. 
arboreus is restricted only to Honshu. We thus conducted 
molecular phylogenetic analyses of the two Rhacophorus 
species from Japan mainland to make their evolutionary 
histories clear. 


2. Materials and Methods 


For the phylogenetic analysis between species to know 
the position of Japanese members among Rhacophorus, 
we employed a total of 17 species, including five species 
from Japan and 12 species from Myanmar, Thailand, 
Malaysia, Vietnam, and China. We determined sequences 
of 12S rRNA, tRNA“, 16S rRNA, and Cytochrome C 
oxidase subunit 1 (CO1) regions of mitochondrial DNA 
for each one sample of each taxon. For outgroup, we 
used each one sample of Polypedates megacephalus in 
the subfamily Rhacophorinae, Buergeria buergeri in the 
subfamily Buergerinae, and Mantella madagascariensis 
of Mantellidae. 

For the intraspecific phylogenetic analyses we 
determined sequences of CO1 for 229 samples of 
R. schlegelii from 156 localities and 95 samples of 
R. arboreus from 70 localities, encompassing whole 
ranges of their distributions (Table 1, Figure 1). For 
the outgroup in the analyses of R. schlegelii, R. viridis, 
R. amamiensis, and R. arboreus were chosen, and for 
those of R. arboreus, R. viridis, R. amamiensis, and R. 
schlegelii were chosen. Additionally, sequences of 12S 
rRNA, tRNA“, and 16S rRNA were determined for R. 
schlegelii, using representative haplotypes obtained in the 
intraspecific analyses of CO1, and utilized for estimation 
of divergence times. 
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We obtained total genomic DNA from the liver or 
muscle samples frozen or preserved in 99% ethanol 
using standard phenol-chloroform extraction procedures 
(Hillis et al., 1996). We conducted amplifications of 
regions shown above by the polymerase chain reaction 
(PCR). The reaction conditions were initial heating 
at 94°C for 4 min; 33 cycles of 94°C (30 s), 53-55°C 
(30 s), and 72°C (2.5 min); and a final extension at 
72°C for 7 min. The amplified DNA fragments were 
purified using polyethylene glycol (PEG, 13%). Cycle- 
sequencing reactions were done using an ABI PRISM 
Big Dye Terminator v3.1 Cycle Sequencing Kit (Applied 
Biosystems) using primers shown in Table 2. We purified 
the sequencing reaction products by ethanol precipitation 
following the manufacturer’s protocol and then ran them 
on an ABI PRISM 3130 Genetic Analyzer (Applied 
Biosystems). 

We aligned sequences using the ClastalW (Thompson 
et al., 1994) or MAFFT version 7 (Katoh and Standley, 
2013), and checked the initial alignments by eye and 
adjusted slightly. We constructed phylogenetic trees by 
maximum likelihood (ML) and Bayesian (BI) methods. 
The optimum substitution models for each partition 
were selected by Kakusan4 (Tanabe, 2011), based on the 
Akaike information criterion. We conducted ML analysis 
using TREEFINDER ver. Mar. 2011 (Jobb, 2011), and 
Bayesian analysis using MRBAYES v3.1.2 (Ronquist 
and Huelsenbeck, 2003) with the models derived from 
Kakusan4 based on the AIC using four simultaneous 
Metropolis-coupled Monte Carlo Markov chains for 
7 000 000 generations and sampling a tree every 100 
generations. Convergence of the log-likelihood scores 
was checked using TRACER ver. 1.5 (Rambaut, 2009), 
and discarded the first 700 000 generations as burn-in. 
We tested the robustness of the ML tree using bootstrap 
analyses (Felsenstein, 1985) with 1000 replicates. We 
regarded tree topologies with bootstrap values (BS) > 
70% as sufficiently resolved (Huelsenbeck and Hillis, 
1993) and posterior probabilities (BPP) for each branch 
in the Bayesian tree > 95% as significant support 
(Huelasenbeck et al., 2001; Leaché and Reeder, 2002). 
We calculated genetic distances of CO1 based on mean 
uncorrected p-distances for pairwise combinations of 
clades and subclades using MEGA, version 5 (Tamura et 
al., 2011). 

For a total 46 samples including 43 samples of the 
two species whose sequences of 12S rRNA, tRNA”, 16S 
rRNA, and CO1 were determined, and three outgroup 
samples, we estimated the divergence times using a 
Bayesian relaxed molecular clock (Drummond et al., 
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Table 1 Samples used for mtDNA analysis in this study with the information of voucher and collection locality. T: Mr. Tanabe's private 
collection. UN: Unnumbered. 


GenBank accession 


. ; 12S rRNA, 
Species Sample No. Locality Voucher (KUHE) Haplotype No. col (RNA, 
16S rRNA 
Rhacophorus schlegelii 1 Mutsu-shi, Aomori Pref. 18149 A-12 LC386398 LC386540 
do. 2 Namioka-machi, Aomori Pref. 29100 A-13 LC386399 - 
do. 3 Outsuchi-cho, Iwate Pref. 29250 A-10 LC386396 LC386541 
do. 4 Oushu-shi, Iwate Pref. 32880 A-10 - - 
do. 5 Ugohonjo-shi, Akita Pref. 31318 A-12 - - 
do. 6 Kesennuma-shi, Miyagi Pref. 45606 A-10 - - 
do. 7 Tomiya-machi, Miyagi Pref. 45609 A-11 LC386397 - 
do. 8 Sendai-shi, Miyagi Pref. 20130827-1 A-11 - - 
do. 9-a Kaminoyama-shi, Yamagata Pref. 20130721-1 A-11 LC386485 - 
do. 9-b Kaminoyama-shi, Yamagata Pref. 20130721-6 A-15 LC386401 - 
do. 10 Takahata-machi, Yamagata Pref. 20130721-33 A-15 - - 
do. 11 Aga-machi, Niigata Pref. 39556 A-15 LC386489 - 
do. 12 Kitakata-shi, Fukushima Pref UN A-15 - - 
do. 13 Aizuwakamatsu-shi, Fukushima Pref. 43338 A-14 LC386400 - 
do. 14 Uonuma-shi, Niigata Pref. 36118 A-5 LC386391 - 
do. ie) Twaki-shi, Fukushima Pref. 21558 A-16 LC386402 - 
do. 16 Utsunomiya-shi, Tochigi Pref. 45604 A-2 LC386388 - 
do. Ly Kamuma-shi, Tochigi Pref. UN A-6 LC386392 - 
do. 18 Midori-shi, Gunma Pref. 31324 A-6 - - 
do. 19 Tsukuba-shi, Ibaraki Pref. 46341 A-4 LC386390 - 
do. 20 Matsudo-shi, Chiba Pref. 13090 A-9 LC386395 LC386542 
do. 21 Ichihara-shi, Chiba Pref. 13095 A-9 - - 
do. 22 Shinano-machi, Nagano Pref. 30680 A-6 - - 
do. 23-a Nagano-shi, Nagano Pref. 46548 A-6 - - 
do. 23-b Nagano-shi, Nagano Pref. 46549 A-6 - - 
do. 23-c Nagano-shi, Nagano Pref. 46550 A-7 LC386393 - 
do. 23-d Nagano-shi, Nagano Pref. 46551 A-7 - - 
do. 23-e Nagano-shi, Nagano Pref. 46552 A-6 - - 
do. 24 Suzaka-shi, Nagano Pref. 40432 A-6 - - 
do. 25 Higashiazuma-machi, Gunma Pref. 30592 A-8 LC386394 - 
do. 26 Ueda-shi, Nagano Pref. 33232 A-3 LC386389 - 
do. 27 Sakuho-machi, Nagano Pref. 46405 A-6 - - 
do. 28 Nagatoro-machi, Saitama Pref. 46342 A-6 - - 
do. 29 Yokoze-machi, Saitama Pref. 20120731-1 A-9 - - 
do. 30-a Atsugi-shi, Kanagawa Pref. 45560 A-6 - - 
do. 30-b Atsugi-shi, Kanagawa Pref. 45614 A-6 - LC386543 
do. 31 Izu-shi, Shizuoka Pref. 20130605-1 A-1 LC386387 LC386544 
do. 32 Itoigawa-shi, Niigata Pref. 30697 D-1-21 LC386448 - 
do. 33-a Kiso-mura, Nagano Pref. 47976 D-1-22 LC386449 - 
do. 33-b Kiso-mura, Nagano Pref. 47977 D-1-22 - - 
do. 34-a Masumigaoka, Ina-shi, Nagano Pref. 47396 D-1-21 - - 
do. 34-b Masumigaoka, Ina-shi, Nagano Pref. 47397 D-1-21 - - 
do. 35-a Takato-machi, Ina-shi, Nagano Pref. 47382 D-1-21 - - 
do. 35-b Takato-machi, Ina-shi, Nagano Pref. 47383 D-1-20 LC386447 - 
do. 36-a Yamanashi-shi, Yamanashi Pref. 20130614-15 D-1-23 LC386450 - 
do. 36-b Yamanashi-shi, Yamanashi Pref. 20130614-28 D-1-23 - - 
do. 37 Fuefuki-shi, Yamanashi Pref. 20130607-1 D-1-21 - - 
do. 38-a Oshino-mura, Yamanashi Pref. 20130614-1 D-1-21 - - 
do. 38-b Oshino-mura, Yamanashi Pref. 20130614-12 D-1-21 - - 
do. 39-a Fujikawa-cho, Yamanashi Pref. 45611 D-1-21 - LC386545 
do. 39-b Fujikawa-cho, Yamanashi Pref. 45612 D-1-21 - 
do. 40-a lida-shi, Nagano Pref. 47402 D-1-21 - - 
do. 40-b lida-shi, Nagano Pref. 47403 D-1-10 LC386438 - 
do. 41 Mori-machi, Shizuoka Pref. 39975 D-1-17 LC386445 - 


do. 42 Ena-shi, Gifu Pref. 30378 D-1-9 LC386437 - 
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Continued Table 1 
GenBank accession 

Species Sample No. Locality Voucher (KUHE) Haplotype No. Col 12S rRNA, 
tRNA“, 16S rRNA 

do. 43-a Shinshiro-shi, Aichi Pref. 26231 D-1-18 LC386446 - 

do. 43-b Shinshiro-shi, Aichi Pref. 26232 D-1-18 - - 

do. 43-c Shinshiro-shi, Aichi Pref. 46067 D-1-18 - LC386546 

do. 44 Seto-shi, Aichi Pref. 48009 D-1-13 LC386441 - 

do. 45 Okazaki-shi, Aichi Pref. 46065 D-1-9 - - 

do. 46 Minamichita-cho, Aichi Pref. 48008 D-1-14 LC386442 - 

do. 47 Nanao-shi, Ishikawa Pref. 28815 D-2a-15 LC386465 - 

do. 48 Kurobe-shi, Toyama Pref. 30700 D-2a-15 - - 

do. 49-a Toga-mura, Nanto-shi, Toyama Pref. 29079 D-2a-15 - LC386547 

do. 49-b Toga-mura, Nanto-shi, Toyama Pref. 29082 D-2a-15 - 

do. 50 Kamitaita-mura, Nanto-shi, Toyama Pref. 13138 D-1-12 LC386440 - 

do. 51 Hakusan-shi, Ishikawa Pref. 33016 D-1-7 LC386435 - 

do. 52 Kamitakara-mura, Takayama-shi, Gifu Pref. 43013 D-1-12 - - 

do. 53-a Shokawa-mura, Takayama-shi, Gifu Pref. 29092 D-1-12 - - 

do. 53-b Shokawa-mura, Takayama-shi, Gifu Pref. 20130606-43 D-1-12 - - 

do. 54 Gujo-shi, Gifu Pref. 16089 D-1-11 LC386439 - 

do. 55 Mino-shi, Gifu Pref. 20130606-38 D-1-12 - - 

do. 56-a Ibigawa-cho, Gifu Pref. 25603 D-1-9 - LC386548 

do. 56-b Ibigawa-cho, Gifu Pref. 25604 D-1-8 LC386436 - 

do. 57 Seki-shi, Aichi Pref. 20131004-2 D-1-15 LC386443 - 

do. 58-a Eiheiji-cho, Fukui Pref. 45830 D-1-21 - - 

do. 58-b Eiheiji-cho, Fukui Pref. 45831 D-1-6 LC386434 - 

do. 58-c Eiheiji-cho, Fukui Pref. 45832 D-1-3 LC386431 - 

do. 59 Fukui-shi, Fukui Pref. 35508 D-1-1 LC386430 - 

do. 60 Echizen-shi, Fukui Pref. 29379 D-1-1 - - 

do. 61 Minamiechizen-cho, Fukui Pref. 26556 D-1-3 - - 

do. 62-a Tsuruga-shi, Fukui Pref. 45838 D-1-5 LC386433 - 

do. 62-b Tsuruga-shi, Fukui Pref. 45839 D-1-5 - - 

do. 62-c Tsuruga-shi, Fukui Pref. 45840 D-1-5 - - 

do. 62-d Tsuruga-shi, Fukui Pref. 45841 D-1-5 - - 

do. 63-a Wakasa-cho, Fukui Pref. 45842 D-1-5 - - 

do. 63-b Wakasa-cho, Fukui Pref. 45843 D-1-5 - - 

do. 63-c Wakasa-cho, Fukui Pref. 45844 D-1-4 LC386432 - 

do. 64 Nagahama-shi, Shiga Pref. 45919 D-1-8 - - 

do. 65-a Taga-cho, Shiga Pref. 14801 D-2a-9 LC386459 - 

do. 65-b Taga-cho, Shiga Pref. 14808 D-1-8 - - 

do. 66-a Takashima-shi, Shiga Pref. 43967 D-2a-2 LC386452 - 

do. 66-b Takashima-shi, Shiga Pref. 45771 D-2a-2 - - 

do. 66-c Takashima-shi, Shiga Pref. 45772 D-2a-2 - - 

do. 67 Otsu-shi, Shiga Pref. 30206 D-2a-5 LC386455 - 

do. 68 Kyoto-shi, Kyoto Pref. 45915 D-2a-5 - - 

do. 69-a Hino-cho, Shiga Pref. 30582 D-2a-10 LC386460 - 

do. 69-b Hino-cho, Shiga Pref. 30583 D-2a-7 LC386457 - 

do. 69-c Hino-cho, Shiga Pref. 31276 D-2a-6 LC386456 - 

do. 69-d Hino-cho, Shiga Pref. 31277 D-2a-6 - - 

do. 70-a Koka-shi, Shiga Pref. 34718 D-2a-8 LC386458 - 

do. 70-b Koka-shi, Shiga Pref. 34719 D-2a-8 - - 

do. 71 Kizugawa-shi, Kyoto Pref. 47975 D-2a-2 - - 

do. 72 Takatsuki-shi, Osaka Pref. 21817 D-2a-11 LC386461 - 

do. 73 Ikoma-shi, Nara Pref. 33262 D-2a-11 - - 

do. 74 Nara-shi, Nara Pref. 45833 D-2a-2 - - 

do. 75-a Nabari-shi, Mie Pref. 31974 D-2a-8 - - 

do. 75-b Nabari-shi, Mie Pref. 31988 D-2a-2 - - 

do. 76 Toshijima-island, Toba-shi, Mie Pref. 42899 D-1-16 LC386444 - 

do. 77 Toba-shi, Mie Pref. 42942 D-2a-14 LC386464 - 

do. 78 Odai-cho, Mie Pref. UN D-2a-12 LC386462 - 

do. 79 Kawakami-mura, Nara Pref. 32729 D-2a-12 - - 
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Continued Table 1 


GenBank accession 
. : 12S rRNA, 
Species Sample No. Locality Voucher (KUHE) Haplotype No. Col tRNA" 
16S rRNA 
do. 80 Hashimoto-shi, Wakayama Pref. 28875 D-2a-12 - - 
do. 8l-a Tomogashima-island, Wakayama-shi, Wakayama Pref. 35994 D-2a-12 - - 
do. 81-b Tomogashima-island, Wakayama-shi, Wakayama Pref. 35996 D-2a-12 - - 
do. 82 Hidaka-cho, Wakayama Pref. 30379 D-2a-12 - - 
do. 83 Totsukawa-mura, Nara Pref. 26251 D-2a-12 - - 
do. 84 Kumano-shi, Mie Pref. 27148 D-2a-13 LC386463 - 
do. 85 Kushimoto-cho, Wakayama Pref. 30384 D-2a-12 - - 
do. 86 Miyazu-shi, Kyoto Pref. 16050 D-2a-2 - LC386550 
do. 87 Kyotanba-cho, Kyoto Pref. 36295 D-2a-1 LC386451 - 
do. 88 Kyotango-shi, Kyoto Pref. 24512 D-2a-4 LC386454 - 
do. 89-a Toyooka-shi, Hyogo Pref. 47597 D-2a-2 - - 
do. 89-b Toyooka-shi, Hyogo Pref. 47598 D-2a-2 - - 
do. 89-c Toyooka-shi, Hyogo Pref. 47599 D-2a-2 - - 
do. 90 Asago-shi, Hyogo Pref. TN2968 D-2a-3 LC386453 - 
do. 91-a Tanba-shi, Hyogo Pref. 47508 B-2 LC386404 - 
do. 91-b Tanba-shi, Hyogo Pref. 47509 B-1 LC386403 - 
do. 91-c Tanba-shi, Hyogo Pref. 47510 B-1 - - 
do. 92-a Awaji-shi, Hyogo Pref. 30052 B-5 LC386407 - 
do. 92-b Awaji-shi, Hyogo Pref. 46070 B-5 - - 
do. 93-a Himeji-shi, Hyogo Pref. 47502 B-10 LC386412 - 
do. 93-b Himeji-shi, Hyogo Pref. 47503 B-4 LC386406 - 
do. 93-c Himeji-shi, Hyogo Pref. 47505 B-7 LC386409 - 
do. 94 Sayo-cho, Hyogo Pref. 36459 B-9 LC386411 - 
do. 95 Tottori-shi, Tottori Pref. 30543 B-6 LC386408 - 
do. 96 Kurayoshi-shi, Tottori Pref. 32799 B-3 LC386405 - 
do. 97 Kofu-cho, Tottori Pref. 45981 B-3 - - 
do. 98-a Hino-cho, Tottori Pref. 45532 B-3 - - 
do. 98-b Hino-cho, Tottori Pref. 45533 B-3 - - 
do. 98-c Hino-cho, Tottori Pref. 45534 B-3 - - 
do. 99-a Nichinan-cho, Tottori Pref. 32790 B-8 LC386410 - 
do. 99-b Nichinan-cho, Tottori Pref. 32791 B-12 LC386414 - 
do. 100 Maniwa-shi, Okayama Pref. 45918 B-3 - - 
do. 101 Tsuyama-shi, Okayama Pref. 45979 B-8 - - 
do. 102 Mimasaka-shi, Okayama Pref. 45978 B-9 - - 
do. 103 Okayama-shi, Okayama Pref. 45529 B-17 LC386419 - 
do. 104 Takahashi-shi, Okayama Pref. 45982 B-12 - - 
do. 105 Niimi-shi, Okayama Pref. 45531 B-3 - LC369670 
do. 106 Jinsekikogen-cho, Hiroshima Pref. 45530 B-11 LC386413 - 
do. 107 Miyoshi-shi, Hiroshima Pref. 45601 B-15 LC386417 - 
do. 108 Mihara-shi, Hiroshima Pref. 30063 B-15 - - 
do. 109 Higashihiroshima-shi, Hiroshima Pref. 45603 B-15 - - 
do. 110 Akiota-cho, Hiroshima Pref. 45600 B-16 LC386418 - 
do. 111 Hiroshima-shi, Hiroshima Pref. 28985 B-14 LC386416 LC386551 
do. 112 Okawa-cho, Kagawa Pref. 30058 B-19 LC386421 - 
do. 113-a Mannou-cho, Kagawa Pref. 30059 B-18 LC386420 - 
do. 113-b Mannou-cho, Kagawa Pref. 30060 B-18 - - 
do. 114 Kamiyama-cho, Tokushima Pref. 22307 B-19 - - 
do. 115-a Aki-shi, Kochi Pref. 47303 B-21 LC386423 - 
do. 115-b Aki-shi, Kochi Pref. 47304 B-19 - - 
do. 115-c Aki-shi, Kochi Pref. 47305 B-19 - - 
do. 116 Kitagawa-mura, Kochi Pref. 29704 B-22 LC386424 LC386552 
do. 117-a Otoyo-cho, Kochi Pref. 47281 B-19 - - 
do. 117-b Otoyo-cho, Kochi Pref. 20130510-1 B-19 - - 
do. 118 Nankoku-shi, Kochi Pref. 47295 B-19 - - 
do. 119 Tosa-shi, Kochi Pref. 47298 B-21 - - 


do. 120-a Nakatosa-cho, Kochi Pref. 47309 B-21 - - 
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Continued Table 1 


GenBank accession 


Species Sample No. Locality Voucher (KUHE) Haplotype No. Col 12S rRNA, 
tRNA“, 16S rRNA 

do. 120-b Nakatosa-cho, Kochi Pref. 47310 B-21 - - 

do. 120-c Nakatosa-cho, Kochi Pref. 47311 B-21 - - 

do. 121 Imabari-shi, Ehime Pref. 47288 B-13 LC386415 - 

do. 122 Kumakogen-cho, Ehime Pref. 47292 B-13 - - 

do. 123 Yusuhara-Cho, Kochi Pref. 47294 B-21 - - 

do. 124 Uwajima-shi, Ehime Pref. 45706 B-19 - - 

do. 125 Tosashimizu-shi, Kochi Pref. 30651 B-20 LC386422 - 

do. 126 Okinoshima-cho, Shimane Pref. 11735 D-2b-1 LC386466 LC386553 
do. 127 Yasugi-shi, Shimane Pref. 20120429-1 D-2b-2 LC386467 - 

do. 128 Unnan-shi, Shimane Pref. 30318 D-2b-3 LC386468 - 

do. 129 Izumo-shi, Shimane Pref. 38804 D-2b-4 LC386469 LC386554 
do. 130-a Oda-shi, Shimane Pref. 45589 D-2b-8 LC386473 - 

do. 130-b Oda-shi, Shimane Pref. 45590 D-2b-8 - - 

do. 131 Kitahiroshima-cho, Hiroshima Pref. 45599 D-2b-8 - - 

do. 132-a Hamada-shi, Shimane Pref. 45594 D-2b-8 - - 

do. 132-b Hamada-shi, Shimane Pref. 45595 D-2b-8 - - 

do. 132-c Hamada-shi, Shimane Pref. 45596 D-2b-8 - - 

do. 133 Masuda-shi, Shimane Pref. 30192 D-2b-13 LC386478 - 

do. 134 Tsuwano-cho, Shimane Pref. 29353 D-2b-14 LC386479 LC386555 
do. 135 Hagi-shi, Yamaguchi Pref. 46003 D-2b-14 - - 

do. 136-a Kuga-machi, Yamaguchi Pref. 30308 D-2b-9 LC386474 - 

do. 136-b Kuga-machi, Yamaguchi Pref. 30309 D-2b-9 - - 

do. 136-c Kuga-machi, Yamaguchi Pref. 30310 D-2b-9 - - 

do. 137 Kitakyushu-shi, Fukuoka Pref. 45903 D-2b-18 LC386483 - 

do. 138 Wakamiya-shi, Fukuoka Pref. 30306 D-2b-17 LC386482 - 

do. 139 Tsukushino-shi, Fukuoka Pref. 20120519-1 D-2b-16 LC386481 - 

do. 140 Saga-shi, Saga Pref. T3195 D-2b-19 LC386484 - 

do. 141-a Tkishima-island, Iki-shi, Nagasaki Pref. 47338 D-2b-15 LC386480 - 

do. 141-b Ikishima-island, Iki-shi, Nagasaki Pref. 47339 D-2b-15 - - 

do. 142 Kawatana-cho, Nagasaki Pref. 17987 D-2b-17 - - 

do. 143-a Nagasaki-shi, Nagasaki Pref. 47227 D-2b-6 LC386471 - 

do. 143-b Nagasaki-shi, Nagasaki Pref. 47228 D-2b-6 - - 

do. 143-c Nagasaki-shi, Nagasaki Pref. 47229 D-2b-6 - - 

do. 144-a Fukuejima-island, Goto-shi, Nagasaki Pref. 31954 C-3 LC386427 - 

do. 144-b Fukuejima-island, Goto-shi, Nagasaki Pref. 31955 C-4 LC386428 

do. 144-c Fukuejima-island, Goto-shi, Nagasaki Pref. 38746 C-5 LC386429 LC386556 
do. 144-d Fukuejima-island, Goto-shi, Nagasaki Pref. 38748 C-4 - - 

do. 144-e Fukuejima-island, Goto-shi, Nagasaki Pref. 45358 C-4 - - 

do. 145 Oita-shi, Oita Pref. 30314 D-2b-17 - - 

do. 146 Arao-shi, Kumamoto Pref. 46069 D-2b-17 - - 

do. 47-a Nishihara-mura, Kumamoto Pref. 20120713-1 C-1 LC386425 LC386557 
do. 47-b Nishihara-mura, Kumamoto Pref. 20120713-7 C-1 - - 

do. 147-c Nishihara-mura, Kumamoto Pref. 20120713-18 C-1 - - 

do. 47-d Nishihara-mura, Kumamoto Pref. 20130606-1 C-2 LC386426 - 

do. 47-e Nishihara-mura, Kumamoto Pref. 20130606-5 C-1 - - 

do. 147-f Nishihara-mura, Kumamoto Pref. 20130606-9 C-1 - - 

do. 47-g Nishihara-mura, Kumamoto Pref. 20130606-11 D-2b-11 LC386476 - 

do. 47-h Nishihara-mura, Kumamoto Pref. 20130606-20 C-1 - - 

do. 48-a Gokase-cho, Miyazaki Pref. 12998 D-2b-17 - LC386558 
do. 148-b Gokase-cho, Miyazaki Pref. 47365 D-2b-17 - - 

do. 148-c Gokase-cho, Miyazaki Pref. 47366 D-2b-17 - - 

do. 149 Amakusa-shi, Kumamoto Pref. 20120426-1 D-2b-14 - - 

do. 150 Satsuma-cho, Kagoshima Pref. 30305 D-2b-7 LC386472 - 

do. 151 Miyazaki-shi, Miyazaki Pref. 30299 D-2b-10 LC386475 - 

do. 152 Miyakonojo-shi, Miyazaki Pref. 30919 D-2b-5 LC386470 - 

do. 153-a Kushima-shi, Miyazaki Pref. 47600 D-2b-5 - - 


do. 153-b Kushima-shi, Miyazaki Pref. 47601 D-2b-5 - - 
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do. 154 Minamiosumi-cho, Kagoshima Pref. 44625 D-2b-12 LC386477 - 
do. 155-a Minamikyushu-shi, Kagoshima Pref. 46988 D-2b-12 - - 
do. 155-b Minamikyushu-shi, Kagoshima Pref. 46989 D-2b-12 - - 
do. 156 Ibusuki-shi, Kagoshima Pref. 30311 D-2b-12 - - 
R. arboreus 1 Mutsu-shi, Aomori Pref. 18091 E-9 LC386494 - 
do. 2 Aomori-shi, Aomori Pref. 29097 E-9 - - 
do. 3 Hirosaki-shi, Aomori Pref. 46309 E-9 - LC386559 
do. 4 Tono-shi, Iwate Pref. 31304 E-9 - - 
do. 5 Oshu-shi, Iwate Pref. 40004 E-10 LC386495 - 
do. 6 Misato-cho, Akita Pref. 40028 E-9 - - 
do. T Yurihonjo-shi, Akita Pref. 40027 E-8 LC386493 - 
do. 8 Oguni-machi, Yamagata Pref. 18464 E-9 - - 
do. 9-a Sado-island, Sado-shi, Niigata Pref. 30809 E-1 LC386485 LC386560 
do. 9-b Sado-island, Sado-shi, Niigata Pref. 30811 E-1 - - 
do. 9-c Sado-island, Sado-shi, Niigata Pref. 30812 E-1 - - 
do. 10 Tadami-cho, Fukushima Pref. 48128 E-7 LC386492 LC386561 
do. 11 Nasushiobara-shi, Tochigi Pref. UN E-4 LC386489 - 
do. 12 Midori-shi, Gunma Pref. 31323 E-4 - - 
do. 13 Hachioji-Shi, Tokyo UN E-4 - - 
do. 14 Joetsu-shi, Niigata Pref. 40333 E-6 LC386491 - 
do. 15-a Itoigawa-shi, Niigata Pref. 30691 E-4 - - 
do. 15-b Itoigawa-shi, Niigata Pref. 30692 E-2 LC386486 - 
do. 15-c Itoigawa-shi, Niigata Pref. 30695 E-2 LC386487 - 
do. 16-c Kurobe-shi, Toyama Pref. 45022 E-4 - - 
do. 16-a Kurobe-shi, Toyama Pref. 30698 E-5 LC386490 - 
do. 16-b Kurobe-shi, Toyama Pref. 30699 E-4 - - 
do. 17-a Shinano-machi, Nagano Pref. 30678 E-3 LC386488 - 
do. 17-b Shinano-machi, Nagano Pref. 46498 E-5 - - 
do. 18 Nagano-shi, Nagano Pref. UN E-4 - - 
do. 19 Azumino-shi, Nagano Pref. UN F-1-15 LC386510 - 
do. 20 Nirasaki-shi, Yamanashi Pref. UN F-1-3 LC386498 - 
do. 21 Minami-alps-shi, Yamanashi Pref. UN F-1-3 - - 
do. 22 Fujikawaguchiko-machi, Yamanaashi Pref. UN F-1-1 LC386496 - 
do. 23 Minobu-cho, Yamanashi Pref. UN F-1-3 - - 
do. 24 Kawazu-cho, Shizuoka Pref. 39969 F-1-2 LC386497 - 
do. 25-a Kanagawa Pref. 47530 - LC386519 - 
do. 25-b Kanagawa Pref. 47531 - LC386520 - 
do. 25-c Kanagawa Pref. 47532 - LC386521 - 
do. 26-a Kimitsu-shi, Chiba Pref. 46165 F-1-4 LC386499 - 
do. 26-b Kimitsu-shi, Chiba Pref. 46166 F-1-4 - - 
do. 27-a Izu-cho, Tokyo Pref. 31938 - LC386522 - 
do. 27-b Izu-cho, Tokyo Pref. 31939 - LC386523 - 
do. 28 Fujieda-shi, Shizuoka Pref. 29349 F-1-5 LC386500 - 
do. 29-a lida-shi, Nagano Pref. 47936 F-1-5 - - 
do. 29-b lida-shi, Nagano Pref. 47945 F-1-5 - LC386562 
do. 30 Urugi-mura, Nagano Pref. UN F-1-5 - - 
do. 31 Neba-mura, Nagano Pref. UN F-1-5 - - 
do. 32 Shitara-cho, Aichi Pref. 28890 F-1-13 LC386508 - 
do. 33 Suzu-shi, Ishikawa Pref. 31302 F-1-14 LC386509 - 
do. 34 Nanto-shi, Toyama Pref. 29056 F-1-12 LC386507 - 
do. 35-a Kanazawa-shi, Ishikawa Pref. UN F-1-15 - - 
do. 35-b Kanazawa-shi, Ishikawa Pref. UN F-1-15 - - 
do. 36 Syokawa-mura, Takayama-shi, Gifu Pref. 29083 F-1-12 - - 
do. 37 Eiheiji-cho, Fukui Pref. 29347 F-1-12 - - 
do. 38 Ibigawa-cho, Gifu Pref. 25605 F-1-13 - - 
do. 39 Nagahama-shi, Shiga Pref. 41473 F-1-12 - - 
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do. 40 Takashima-shi, Shiga Pref. 29291 F-1-12 - - 
do. 41 Maibara-shi, Shiga Pref. 46419 F-1-12 - - 
do. 42 Kameyama-shi, Mie Pref. UN F-1-11 LC386506 - 
do. 43-a Odai-cho, Mie Pref. 30826 F-1-12 - - 
do. 43-b Odai-cho, Mie Pref. 30827 F-1-16 LC386511 - 
do. 44 Ujitawara-cho, Kyoto Pref. 45068 F-1-12 - - 
do. 45 Kyoto-shi, Kyoto Pref. 45069 F-1-12 - LC386563 
do. 46 Kyotanba-cho, Kyoto Pref. 46466 F-1-12 - - 
do. 47 Kameoka-shi, Kyoto Pref. T3479 F-1-8 LC386503 - 
do. 48 Takatsuki-shi, Osaka Pref. 21813 F-1-12 - - 
do. 49 Miyazu-shi, Kyoto Pref. 16049 F-1-12 - - 
do. 50 Fukuchiyama-shi, Kyoto Pref. 46538 F-1-9 LC386504 - 
do. 51 Kyotango-shi, Kyoto Pref. T3456 F-1-10 LC386505 - 
do. 52-a Toyooka-shi, Hyogo Pref. 30584 F-1-8 - - 
do. 52-b Toyooka-shi, Hyogo Pref. T3055 F-1-9 - - 
do. 53-a Asago-shi, Hyogo Pref. 30532 F-1-9 - - 
do. 53-b Asago-shi, Hyogo Pref. 30533 F-1-8 - - 
do. 53-c Asago-shi, Hyogo Pref. 36585 F-1-8 - - 
do. 54 Sayo-cho, Hyogo Pref. 30537 F-1-7 LC386502 - 
do. 55 Tottori-shi, Tottori Pref. 30542 F-1-8 - - 
do. 56-a Kagamino-cho, Okayama Pref. 30538 F-1-6 LC386501 - 
do. 56-b Kagamino-cho, Okayama Pref. 32820 F-1-8 - - 
do. 56-c Kagamino-cho, Okayama Pref. UN F-1-17 LC386512 - 
do. 57 Daisen-cho, Tottori Pref. 31051 F-2-2 LC386514 - 
do. 58 Niimi-shi, Okayama Pref. 45046 F-2-2 - LC386564 
do. 59 Matsue-shi, Shimane Pref. 31056 F-2-3 LC386515 - 
do. 60-a Saijo-cho, Shobara-shi, Hiroshima Pref. 31057 F-2-2 - - 
do. 60-b Saijo-cho, Shobara-shi, Hiroshima Pref. 31058 F-2-1 LC386513 - 
do. 61 Tojo-cho, Shobara-shi, Hiroshima Pref. 31059 F-2-2 - - 
do. 62 Takano-cho, Shobara-shi, Hiroshima Pref. 31060 F-2-2 - - 
do. 63-a Izumo-shi, Shimane Pref. 31053 F-2-2 - - 
do. 63-b Izumo-shi, Shimane Pref. 31054 F-2-3 - - 
do. 63-c Izumo-shi, Shimane Pref. 31055 F-2-2 - - 
do. 64-a Miyoshi-shi, Hiroshima Pref. 31061 F-2-2 - - 
do. 64-b Miyoshi-shi, Hiroshima Pref. 31062 F-2-4 LC386516 - 
do. 65 Kitahiroshima-cho, Hiroshima Pref. 31064 F-2-5 LC386517 - 
do. 66-a Akiota-cho, Hiroshima Pref. 31063 F-2-3 - - 
do. 66-b Akiota-cho, Hiroshima Pref. 47579 F-2-5 - LC386565 
do. 66-c Akiota-cho, Hiroshima Pref. 47580 F-2-5 - - 
do. 67 Hiroshima-shi, Hiroshima Pref. 28975 F-2-3 - - 
do. 68 Tsuwano-cho, Shimane Pref. 29348 F-2-5 - - 
do. 69 Toguchi-cho, Yamaguchi-shi, Yamaguchi Pref. 28995 F-2-6 LC386518 - 
do. 70 Yamaguchi-shi, Yamaguchi Pref. 28994 F-2-5 - - 
R. cyanopunctatus KotaTinggi, Malaysia 52122 - LC386526 LC386566 
R. dennysi Tam Dao, Vietnam, 55211 - LC386527 LC010575 
R. dugritei Emeishan, Sichuan, China 27701 - LC386528 LC010584 
R. dulitensis Tawau, Malaysia B09087 - LC386529 AB847123 
R. feae Kachin, Myanmar 34511 - LC386530 LC386567 
R. maximus Pilok, Thailand 35130 - LC386531 LC386568 
R. minimus China 70049 - LC386532 LC386569 
R. moltrechti Taipei, Taiwan, China 31070 - LC386533 LC386570 
R. nigropunctatus China 70073 - LC386534 LC386571 
R. norhayati Malaysia UN - LC386535 AB728191 
R. omeimontis Emeishan, Sichuan, China 20060104 - LC386536 LC010595 
R. owstoni Ishigakijima-island, Okinawa Pref. 12764 - LC386537 LC386572 
R. rhodopus Jinghong, China 47794 - LC386538 LC386573 
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A : 12S rRNA, 
Species Sample No. Locality Voucher (KUHE) Haplotype No. col tRNA“, 
16S rRNA 
R. taipeianus Taipei, Taiwan, China 34347 - LC386539 LC386574 
R. amamiensis Amamioshima-island, Kagoshima Pref. 22524 - LC386524 LC386575 
R. viridis Okinawajima-island, Okinawa Pref. 35354 - LC386525 LC386576 
Polypedates megacephalus Guilin, Guangxi, China - NC006408 
Buergeria buergeri Hiroshima Pref., Japan - AB127977 
Mantella madagascariensis Madagascar - AB212225 


2006) calculated in BEAST ver. 1.6.2 (Drummond and 
Rambaut, 2007) using 10 million generations, discarding 
the first 1 million generations as burn-in. Suitable burn- 
in and convergence of parameters were determined using 
TRACER ver. 1.5 (Rambaut, 2009). As calibration points, 
the divergence of Buergerinae and Rhacophorinae as 45.9 
(95% credibility interval [CI] 57.6-36.6) million years 
before present (MYBP: Kim et al., 2007), of Rhacophorus 
and Polypedates as 30.6 (CI 34.7-25.2) MYBP (Li et 
al., 2013), and of Cynops ensicauda from Okinawa and 
Amami Islands as 5.21 (CI 6.82-3.64) MY BP (Tominaga 
et al., 2013) were used. 

We constructed haplotype networks for a set of primers 
(588 bp, Rh_CO1_F1 and Rh_CO1_R1) in CO1 for 
R. schlegelii and R. arboreus, and for each clade and 
subclade, we calculated both haplotype diversity (A) and 
nucleotide diversity (n, based on pairwise differences) 
values and conducted neutrality tests. We used TCS 
ver. 1.21 (Clement et al., 2000) to construct haplotype 
networks. Statistical confidence limit was set at >95%. 
For R. schlegelii, network was constructed for each clade 
and subclade since networks were not constructed in the 
confidence interval of 95%. For calculation of haplotype 
and nucleotide diversities, and conducting neutrality 
tests, ARLEQUIN ver. 3.5.1.3 (Excoffier and Lischer, 
2010) was used. Putative populations of introduced origin 
(Localities 25: Kanagawa Pref. and 27: Izu-Oshima, see 
below in Discussion) were not included in these analyses. 


3. Results 


3.1. Interspecific phylogeny of Rhacophorus Of 3976 
nucleotide sites of 12S rRNA, tRNA”, 16S rRNA, 
and CO1 regions obtained for Rhacophorus including 
Japanese members, 1422 were variable within the 
ingroup. For the ML analysis, J2 (Jobb, 2008)+G was 
selected by Kakusan4 as the optimal model for 12S rRNA 
and tRNA“, J2 (Jobb, 2008)+G for 16S rRNA, and TN93 
(Tamura and Nei, 1993)+G, HKY85 (Hasegawa et al., 
1985)+G, and J3 (Jobb, 2008)+G were selected as the 


optimal models for Ist, 2nd, and 3rd codon positions of 
CO1 gene. For the BI analysis, GTR+G was selected for 
12S rRNA, tRNA”, 16S rRNA, and SYM (Zharkikh, 
1994)+G, HKY85+G, and GTR+G were selected for 
Ist, 2nd, and 3rd codon positions of CO1 gene. The ML 
and Bayesian analyses produced topologies with -In L = 
21 917.505 and 21 918.380, respectively. 

Phylogenetic analyses employing two different 
optimality criteria yielded nearly identical topologies. 
As shown in the ML tree (Figure 2), Japanese and some 
Chinese species formed a clade (ML BS = 99%, BPP 
= 1.00) in which R. owstoni and R. moltrechti formed 
a clade (100%, 1.00), which was sister to the clade 
(70%, 1.00) of other Japanese taxa and Chinese species. 
However, other than the monophyly of R. viridis and R. 
amamiensis, phylogenetic relationships within Japanese 
Rhacophorus were not resolved. 


3.2. Intraspecific phylogeny of R. schlegelii From 
1554 nucleotide sites of complete CO1 region obtained 
for 229 individuals of R. schlegelii, 150 haplotypes 
were recognized, and 420 positions were found variable 
including the outgroup. For the ML analysis, Kakusan4 
selected TN93+G, HKY85, and TN93+G as the optimal 
models for 1st, 2nd, and 3rd codon positions of CO1 
gene, respectively. For the BI analysis, SYM+G, F81 
(Felsenstein, 1981)+G, and GTR+G, were selected 
respectively, for each position 

Phylogenetic analyses yielded identical topologies in 
the two different criteria, and only the ML tree is shown 
(Figure 3A). Rhacophorus schlegelii formed a clade with 
high supports (99%/1.00) in which four major clades, 
each with high support, were recognized (A: 99%/1.00; 
B: 99%/1.00; C: 98%/1.00; D: 93%/1.00), although their 
relationships were unresolved. In addition, three subclades 
were recognized within Clade D (D-1: 99%/1.00; D-2a: 
99%/1.00; D-2b: 100%/1.00). 

In the phylogenetic trees obtained for 12S rRNA, 
tRNA“, 16S rRNA, and CO1, all the four clades obtained 
in the CO1 trees were recognized, and Clade A proved to 
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Figure 1 Map of sampling localities of (A) Rhacophorus schlegelii and (B) R. arboreus. For locality name, see Table 1. 
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Figure 2 Maximum likelihood (ML) tree based on mitochondrial 12S rRNA, tRNA“, 16S rRNA and CO1 genes of the genus Rhacophorus. 
Nodal values indicate bootstrap for ML and Bayesian posterior supports (ML-BS/BSS). For voucher number and DDBJ accession number, 
see Table 1. 


be a sister clade to the clade of the remaining three clades As shown in Figure 4A, range of distribution was 
(87%/0.99), each of which formed a clade (B: 99%/1.00; Tohoku, Kanto, and eastern Chubu regions of Honshu 
C: 99%/1.00; D: 98%/1.00), but their relationships were Island in Clade A, westernmost Kinki and eastern 
unresolved. Chugoku regions of Honshu, and Shikoku Island in 
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Figure 3 Maximum likelihood (ML) tree based on mitochondrial CO1 gene of (A) R. schlegelii and (B) R. arboreus. Nodal values indicate 
bootstrap supports for ML and Bayesian posterior probably (ML-BS/BSS). For locality numbers, see Figure 1 and Table 1. 
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Figure 4 Geographic distribution of each clade and subclade of (A) R. schlegelii and (B) R. arboreus. For sample marks and lines, refer to 


Figure 1. 
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Clade B, and a part of Kyushu Island (localities 144 and 
147) in Clade C. Clade D had wider range than others 
in Chubu, Kinki, western Chugoku and Kyushu regions, 
with Subclade D-1 in Chubu region except for north and 
Toshijima Island (locality 76), Subclade D-2a in northern 
Chubu and Kinki regions, and Subclade D-2b in western 
Chugoku and Kyushu regions. Of these, sympatric 
distribution was found between Clades C and D (at 
locality 147) and Subclades D-1 and D-2a (at locality 65). 
Other clades and subclades occur parapatrically (Figure 
4A). 

Mean (+ ISE, in Yo) of uncorrected p-distance was 
5.09 + 0.51 between Clades A and B, 4.61 + 0.50 between 
Clades A and C, 5.36 + 0.48 between Clades A and D, 
4.70 + 0.47 between Clades B and C, 4.93 + 0.39 between 
Clades B and D, and 4.80 + 0.43 between Clades C and 
D. In Clade D, they were 4.56 + 0.48 between Subclades 
D-1 and D-2a, 3.69 + 0.39 between D-1 and D-2b, and 
3.30 + 0.41 between D-2a and D-2b. 


3.3. Intraspecific phylogeny of R. arboreus In R. 
arboreus, 58 haplotypes were recognized within 95 
individuals. Of 1 554 bp of complete CO1 region, 347 
positions were found variable. For the ML analysis, 
Kakusan4 selected TN93+G, HKY85, and J1 (Jobb, 
2008)+G, respectively, as the optimal models for Ist, 
2nd, and 3rd codon positions of CO1 gene. For the BI 
analysis, SYM+G, F81+G, and GTR+G, were selected 
respectively, for each position. 

The ML and Bayesian analyses produced topologies 
with -ln L = 4 240.711 and 4 247.162, respectively. 
Identical topologies were obtained in ML and BI, and ML 
tree is shown in Figure 3B. Rhacophorus arboreus formed 
a clade with high supports (99%/1.00), within which two 
major clades, Clade E (98%/1.00) and Clade F (86/0.95), 
were recognized, and two Subclades (F-1: 78%/1.00 and 
F-2: 98%/1.00) were recognized in the latter. 

Clade E occurs in Tohoku, Kanto except for southern 
part, and eastern Chubu regions, and Clade F in southern 
Kanto, Chubu, Kinki, and Chugoku regions (Figure 
4B). In Clade F, Subclade F-1 is found in southern 
Kanto, Chubu, and Kinki regions, and Subclade F-2 in 
Chugoku regions. These clades and subclades occurred 
allopatrically. Topotypic samples of synonymized 
subspecies R. schlegelii intermedia (locality 9) were 
included in Subclade F-1 without remarkable genetic 
differentiation among R. arboreus from the surrounding 
area. 


Mean (+ ISE, in Yo) of uncorrected p-distance was 2.40 
+ 0.33 between Clades E and F, and 1.46 + 0.24 between 
F-1 and F-2. 
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3.4. Divergence time Rhacophorus owstoni formed a 
clade with R. moltrechti in Taiwan, China and split from 
the group containing other Japanese Rhacophorus at 
7.68 (9.44-6.05) MYBP. Rhacophorus owstoni and R. 
moltrechti were estimated to have diverged at 2.32 (3.35- 
1.47) MYBP. Japanese members other than R. owstoni 
were estimated to have split from the common ancestor of 
Chinese species including R. omeimontis and R. minimus 
at 7.01 (8.53-5.51) MYBP. Divergence of R. viridis and 
R. amamiensis was estimated to be 3.57 (4.56—2.58) 
MYBP. 

Within R. schlegelii, Clade A and others first diverged 
at 2.40 (3.02-1.79) MYBP, and divergences among 
Clades B, C, and D occurred after 2.04 (2.57—1.52) 
MYBP. Differentiation within each clade and subclade 
was estimated to be 0.51 (0.71-0.32) MYBP in Clade 
A, 0.84 (1.15-0.55) MYBP in CladeB, 0.52 (0.79-0.27) 
MYBP in Clade C, 1.44 (1.80-1.05) MYBP between 
Subclade D-1 and others in Clade D, and 1.20 (1.53-0.87) 
MYBP between Subclades D-2a and D-2b. Differentiation 
within Subclade D-1 was 0.76 (1.02-0.50) MY BP, D-2a 
was 0.66 (0.90-0.42) MY BP, and D-2b was 0.44 (0.64— 
0.25) MYBP. 

Within R. arboreus, divergence of Clades E and F was 
estimated to have occurred at 0.96 (1.30-0.65) MY BP. 
Within each clade or subclade, age of differentiation was 
estimated to be 0.23 (0.35—0.12) MYBP in Clade E, 0.65 
(0.90-0.42) MYBP in Clade F, 0.30 (0.46-0.14) MYBP 
in Subclade F-1, and 0.19 (0.30—0.08) MYBP in Subclade 
F-2. 


3.5. Haplotype network and population genetic 
analysis Haplotype network of R. schlegelii, based on 
588 bp of COI region of each clade and subclade is 
shown in Figure 5A, and haplotype diversity, nucleotide 
diversity, and results of neutrality tests (Fu’s FS: Fu, 
1997) are shown in Table 3. Clade A tended to have lower 
nucleotide diversity than others, and Clades A and C 
tended to have lower haplotype diversity than others. Fu’s 
FS were significantly negative (P < 0.05) in Clades A and 
D, and Subclade D-2b. 

Haplotype network of R. arboreus, based on 588 bp 
of CO1 region of each clade and subclade is shown in 
Figure 5B, and haplotype diversity, nucleotide diversity, 
and result of test of neutrality (Fu’s FS) are shown in 
Table 3. Mean of haplotype diversities (0.8646) was 
slightly smaller and nucleotide diversity (0.0043) was 
much smaller than in R. schlegelii (0.8981 and 0.0124, 
respectively). In the two clades recognized (E and F) and 
Subclade F-1, Fu’s FS were significantly negative, but in 
Subclade F-2, the values were not significantly negative. 
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Table 2 Primers used to amplify mtDNA in this study. 

Target Name Sequence (5'—3') Reference 
12SL1 RYAGTGY TGAAAACGCTAAGA Matsui et al. (2010) 
12SH1 GACACCGTCAAGTCCTTTGGGTTT This study 
12SA-L AAACTGGGATTAGATACCCCACTAT Palumbi et al.(1991) 

12S rRNA and tRNA” H15548 TACCATGTTACGACTTTCCTCTTCT Matsui et al.(2005) 
Lnew TACACACCGCCCGTCACCCTCTT Shimada et al. (2011) 
tVal-H AAGTAGCTCGCTTAGTTTCGG Shimada et al. (2011) 
tVal-L CGTACCTTTTGCATCATGGTC Shimada et al. (2011) 

H2315Po TTCTTGTTACTAATCY YAGCAT Kuraishi et al. (2012) 
L2188 AAAGTGGGCCTAAAAGCAGCCA Matsui et al. (2006) 
16S rRNA Wilkinson_6 CCCTCGTGATGCCGTTGATAC Wilkinson et al. (2002) 
16L1 CTGACCGTGCAAAGGTAGCGTAATCACT Hedges (1994) 
16H1 CTCCGGTCTGAACTCAGATCACGTAGG Hedges (1994) 
Rh CO1 F1 CCTCTGTCTTTGGGGCTACA This study 
Rh CO1 R1 GCAGCTAGTACWGGGAGRGA This study 
col Rh CO1 F2 GGAGCTATTAAYTTTATYACAAC This study 
Rh CO1 R2 TCAGTGGATAAAKCCTGCYATAAT This study 
Rh CO1 F3 GCCAACTCTTCTATTGAYATTGT This study 
Rh CO1 R3 TTAACTTGYAAWTTTRCCTTGACA This study 
Table 3 Genetic diversity indices and result of Fu’s tests of neutrality (Fu, 1997). 
Clade and subclade Number of samples Number of haplotypes Haplotype diversity Nucleotide diversity Fu’s FS Fu’s FS P-value 
A 37 16 0.8919 + 0.0381 0.004402 + 0.002681 -7.92118 0.000 
B 49 22 0.9379 + 0.0173 0.013652 + 0.007159 -3.21333 0.165 
C 12 5 0.7273 0.1133 0.011054 + 0.006334 3.05309 0.925 
R. schlegelii D 131 57 0.9740 + 0.0047 0.029069 + 0.014402 -10.10994 0.049 
D-1 50 23 0.9339 + 0.0202 0.013836 + 0.007245 -3.72552 0.139 
D-2a 39 15 0.8853 + 0.0320 0.009135 + 0.005001 -1.77254 0.275 
D-2b 42 19 0.9361 + 0.0185 0.006232 + 0.003574 -7.89028 0.001 
E 25 10 0.8633 + 0.0440 0.002755 + 0.001878 -4.79453 0.004 
Fa TWA F 65 23 0.9298 + 0.0155 0.007323 + 0.004072 -7.70173 0.009 
F-1 45 17 0.8970 + 0.0295 0.004408 + 0.002671 -8.04489 0.001 
F-2 20 6 0.7684 + 0.0624 0.002900 + 0.001974 -0.71361 0.338 


4. Discussion 


4.1. Phylogenetic relationships of Japanese 
Rhacophorus Li et al. (2013) and Hertwig et al. (2013) 
similarly recognized three clades (I: mainly Bornean 
species; II: flying frogs; III: species occurring mainly 
in East Asia) in the genus Rhacophorus. Our result 
corroborated them and all Japanese Rhacophorus are 
included in the group mainly from East Asia (Clade IID), 
which formed a monophyletic group against Clade I (R. 
cyanopunctatus and R. gauni) and Clade II (R. norhayati 
and R. rhodopus). Within Clade HI of mainly East Asian 
species group, southern species like R. dulitensis, R. 
feae, and R. maximus diverged first and Japanese species 
positioned relatively terminally on the phylogram. From 
this relationship, species in this group are estimated to 
have originated in the south and sometime later invaded 
north to form present distribution. 

It is difficult to elucidate interspecific phylogenetic 


relationships among Rhacophorus, but among Japanese 
taxa, R. owstoni from the Yaeyama group of the southern 
Ryukyus formed a clade with R. moltrechti from Taiwan, 
China and split from the remaining Japanese and Chinese 
species. Except for R. owstoni, Japanese Rhacophorus had 
ancestor common to species like R. taipeianus in Taiwan, 
China and Chinese R. omeimontis, R. dugritei, and R. 
minimus, although their relationships were unresolved. It 
was estimated that the differentiation of these frogs began 
7.01 MYBP, which is much older than the period of last 
disjunction of Ryukyu Archipelago and mainland Japan 
(1.70 MYBP: Kitamura and Kimoto, 2004). From this it 
is estimated that divergence of ancestors of Japan and the 
continent already occurred by unknown factors before 
the isolation of Japanese Islands and Ryukyu Archipelago 
from continental China by oceans. 


4.2. Phylogenetic relationships in R. schlegelii On the 
phylogenetic tree based on multiple genes, R. schlegelii 
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Figure 5 Statistical parsimony network of each clade and subclade of (A) R. schlegelii and (B) R. arboreus. Filled circles indicate missing 
haplotypes. The size of each open circle is proportional to the haplotype frequency. 
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was split into four clades. Among them Clade A was 
the first to diverge from the others, and therefore, R. 
schlegelii was divided into two groups roughly from 
eastern and western regions (Figure 4A). 

Near the boundaries of these clades, Japan Alps and 
Kanto Mountains are situated between Clades A and 
D, and Chugoku Mountains between Clades B and D. 
Current boundary of Clades A and D is positioned in the 
east of Japan Alps, but the boundary is estimated to have 
been more west in the past since the haplotype of Clade 
A occurs isolated on the Izu Peninsula. Clade A and the 
other clades were estimated to have diverged around 2.40 
MYBP. It is estimated that the altitude of Hida Mountains 
(northern Japan Alps) was ca. 2000 m but of Akaishi 
Mountains (southern Japan Alps) was ca. 1000 m in that 
era (Chinzei and Machida, 2001). Therefore, it is not 
likely that Akaishi Mountains at that era were a barrier 
for the migration of R. schlegelii, which mainly inhabits 
from lowland to subalpine regions. However, around 2.70 
MY BP, which nearly corresponds to the age of divergence 
between Clade A and the others (2.40 MYBP), cycles of 
the glacial and the interglacial periods are estimated to 
become evident (Mix et al., 1995), and climatic factors 
like lowering of air temperatures on high mountains 
would have inhibited population gene exchange and 
promoted the differentiation of R. schlegelii. 

Boundary of Clades B and D is roughly situated near 
Chugoku Mountains. On the Japan Sea side, Clade B is 
distributed between Subclades D-2a and D-2b. In the past, 
Subclades D-2a and D-2b are estimated to have occurred 
in the north, while Clade B occurred in the south, of 


Chugoku Mountains. Subsequently, Clade B is estimated 
to have widened the range locally in the north of Chugoku 
Mountains. It is estimated that Subclades D-2a and D-2b 
first differentiated and Clade B later was distributed at 
their boundary, since differentiation within Clade B was 
0.84 MYBP, while divergence of Subclades D-2a and 
D-2b was 1.20 MYBP. It is possible that divergence of 
Subclades D-2a and D-2b was affected by the presence 
of active volcano of Mt. Daisen and cooler climate of 
Japan Sea side, but the volcanic activity of Daisen is 
estimated to have begun 0.90 MYBP (Ota et al., 2004), 
and there seems to have been no significant climatic 
variation. Furthermore, altitudes of Chugoku Mountains 
are at most 800-1000 m a.s.l. (Ota et al., 2010), which is 
hardly considered to be a barrier of distribution, inducing 
differentiations of R. schlegelii. Thus the reasons for 
differentiations of Clades B and D, and of Subclades D-2a 
and D-2b are unknown. 

Clades C and D occurred sympatrically at the 
location 147 (Kumamoto, western central Kyushu). It is 
conceivable that Clade D invaded to the range of Clade C 
in Kyushu and occupied present areas. 

Differentiations within Clade A were less marked than 
in other clades, and haplotype diversity and nucleotide 
diversity were smaller than in Clades B and D (Table 
3). Additionally, results of neutrality tests indicated 
past rapid expansion of distribution range, and the 
haplotype exhibited star-like network (Figure 5A), with 
the commonest haplotype (A-6) occurring in Kanto and 
eastern Chubu districts, indicating that Clade A utilized 


these regions as refugia in the past (Figure 6A). 


Figure 6 Examples of dispersals in (A) R. schlegelii (Clade A [A-6]) and (B) R. arboreus (Subclade F-1 [F-1-12]) that were haplotypes most 


frequently observed. 


No. 2 Masafumi MATSUI et al. 


Degree of genetic differentiation in Clade B, as 
shown by haplotype diversity and nucleotide diversity 
was smaller than Clade D, which includes diversified 
subclades, but was larger than in Clades A and C (Table 
3). Haplotype network contains long branches among 
each haplotype (Figure 3). This indicates that Clade B has 
originated older than Clade A but younger than Clade D, 
and sustained population size without experiencing large- 
scale changes in distribution. 

Clade C from two distant localities in Kyushu had 
lower haplotype diversity than in other clades, but had 
nucleotide diversity higher than Clade A, Subclades D-2a 
and D-2b (Table 3). Divergence time between disjunct 
populations was 0.52 MYBP, and are estimated have 
been distributed widely on Kyushu at least by that period. 
Probably, the clade was substituted by Subclade D-2b, 
which invaded after Clade C decreased the population 
size through volcanic activities in northern Kyushu (2.0— 
0.5 MYBP: Sakamoto, 2009). Occurrence of genetically 
different lineages in Kyushu is known in the newt Cynops 
pyrrhogaster (Tominaga et al., 2013) and the brown frog 
Rana tagoi (Eto et al., 2012), and Clade C might have 
experienced geohistrical events similar to them. 

Clade D showed larger genetic differentiation than 
in other clades, and three fairly divergent subclades 
were recognized. Of these, Subclade D-2b was 
estimated to have experienced rapid range expansion 
based on neutrality test (Table 3). The subclade began 
differentiation at 0.44 MYBP, and is estimated to have 
occurred in westernmost Honshu and northern Kyushu, 
because the dominant haplotypes (D-2b-8, 17) are 
obtained from these areas (Figure 5, Table 1). Then, 
expansion of distribution towards southern Kyushu 
on Subclade D-2b seems to have induced restricted 
distribution of Clade C shown above. 


4.3. Phylogenetic relationships in R. arboreus 
Phylogenetic analyses using CO1 gene revealed the 
presence of eastern Clade E and western Clade F in 
R. arboreus, with further two subclades in Clade F 
(Figures 3B and 4B). These clades and subclades occur 
parapatrically, and the boundary of Clades E and F 
was near the region connecting Japan Alps and Kanto 
Mountains. However, because no differentiation is found 
between many populations with mountains or mountain 
ranges in between, Japan Alps and Kanto Mountains are 
hardly considered to have induced genetic differentiation. 
The degree of genetic differentiation is much smaller than 
in R. schlegelii: mean uncorrected p-distance was 2.4% 
between Clades E and F of R. arboreus, as compared with 
5.4% between Clades A and D of R. schlegelii. Thus, R. 
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arboreus might have experienced events where they lost 
genetic diversities, and retreated to refugia, indicating 
genetic differentiations in R. arboreus to have been 
induced by their retreats to refugia. 

Clade E from eastern Japan had lower degree of 
genetic differentiation, haplotype diversity, and nucleotide 
diversity than western Clade F (Table 3), and was 
suggested to have experienced past rapid expansion of 
distribution range based on the neutrality test. In the range 
of Clade E, haplotypes E-4 in Tohoku and E-9 in northern 
Chubu were found in high frequency and were distributed 
widely. This suggests that the original populations of 
these haplotypes that started their rapid range expansions 
were present at their boundaries, which would have been 
their refugia. 

Within Clade F, the Subclade F-1 showed significant 
negative value in the neutrality tests and suggested past 
rapid range expansion (Table 3). Within Subclade F-1, 
the haplotype F-1—12 was found frequently, and exhibited 
star-like haplotype network (Figure 5B). This haplotype 
occurred widely from Kinki to northern Chubu regions, 
and might utilize this area as their refugia (Figure 6B). 

No evidence of range expansion or reduction was 
obtained in the Subclade F-2, and differentiation in the 
subclade was 0.19 MYBP, in contrast to 0.65 MYBP 
estimated for the differentiation from Subclade F-1. 
The haplotype network exhibits array-like pattern. This 
indicates that the Subclade F-2 was already distributed 
in the Chugoku Mountains in 0.65 MYBP and sustained 
stable population size and distributional range. In this 
way, R. arboreus seems to have experienced to retreat to 
the refugia, and caused genetic differentiations among 
local populations thus formed. It is estimated that the 
populations in central Honshu retreated to northern 
regions, as if avoiding warm climates, but no evidence 
was found for such a hypothesis. Retreats to refugia might 
relate to the limited distribution of R. arboreus to Honshu, 
but we could not get any evidence for that. 

Localities 25 (Kanagawa Pref.) and 27 (Izu-Oshima) 
were isolated from other localities in Clade F. Of these, 
locality 27 is represented by artificially introduced 
populations (Maeda and Matsui, 1999), and locality 
25, with the haplotype same as Kinki region, is also 
considered an alien population. 


4.4. Comparative phylogeography This study clarified 
that both R. schlegelii (Clade A and others) and R. 
arboreus (Clades E and F) are genetically divergent 
between the east and west regions of Japan. Such eastern 
and western genetic divergence is known in several 
wide-ranging Japanese frogs, although the boundaries 
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are variable (e.g., Bufo japonicus [Igawa et al., 2006]; 
Hyla japonica [Dufresnes et al., 2016]; Rana japonica 
[Sumida and Ogata, 1998]). The east-west boundary of 
R. schlegelii positioned in the east of Japan Alps, which 
is common to R. arboreus, except for the distribution 
of Clade F in Chubu and southern Kanto regions. 
Additionally, the boundary of Subclades D-2a and D-2b 
in R. schlegelii positioned in Chugoku region, like the 
boundary of Subclades F-1 and F-2 in R. arboreus. In 
this way, boundaries of distributions between genetic 
groups corresponded well in the two Rhacophorus 
species. However, estimated ages of differentiations 
greatly differed between the two species and their current 
distributions are hardly attributable to identical geological 
events or geographic boundaries. 

Both the two species of Rhacophorus were suggested 
to include genetic populations that experienced 
retreatment to refugia. Among them, typical cases are 
as follows: in R. schlegelii, populations now occurring 
in Tohoku, Kanto, and eastern Chubu regions (Clade A) 
are estimated to have utilized Kanto and eastern Chubu 
regions as their refugia, while in R. arboreus, populations 
occurring Chubu to Kinki regions (Subclade F-1) are 
estimated to have refugia in the northern part of that 
regions. In addition, Subclade D-2b in R. schlegelii and 
Clade E in R. arboreus indicate the trace of retreat to 
refugia. 

Rhacophorus arboreus was found much less 
diversified genetically than R. schlegelii, although it is 
more variable than the latter species in morphology and 
breeding habits (Nakamura and Uéno, 1963; Matsui and 
Maeda, 2018). At present, Clade A of R. schlegelii and 
Clade E of R. arboreus occur sympatrically, but because 
R. schlegelii is estimated to have retreated to Kanto and 
eastern Chubu regions, wide areas of eastern Japan is 
estimated to be occupied only by R. arboreus. In western 
Japan, R. arboreus and R. schlegelii are estimated to 
have coexisted long because R. arboreus retreated to 
north of Kinki and Chubu regions. In the regions where 
the two species co-occurred, they interfered each other, 
which might have resulted in enhancement of the habit 
of breeding on trees and possession of body markings in 
the western populations of R. arboreus. In order to test 
this hypothesis, knowledge of morphology and breeding 
habits, both now considered not variable, of R. schlegelii 
should be accumulated and interspecific interference at 
the sympatric zone should be studied ecologically. 
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